function y = Cd_a(p, x, z)

  if x(12) < -p.b  % ground effect
      c = [0 558.19 -212.69 17.655 2.0913 -0.0559 0.0227];
      e = [z.^ 6, z.^ 5, z.^ 4, z.^ 3, z.^ 2, z.^ 1, z.^ 0]';
  else
      c = [0 0 0 -0.6809 2.141 -0.0606 0.0263];
      e = [z.^ 6, z.^ 5, z.^ 4, z.^ 3, z.^ 2, z.^ 1, z.^ 0]';
  end
  
  y = c * e;
  
end